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ABSTRACT 

We present time-resolved optical spectroscopy of the old nova RRPic. Two emission 
lines (Ha and He I) are present in the observed part of the spectrum and both show 
strong variability. Ha has been used for Doppler tomography in order to map the 
emission distribution in this system for the first time. The resulting map shows the 
emission from the disc as well as two additional emission sources on the leading and 
trailing side of the disc. Furthermore we find evidence for the presence of either a disc- 
overflow or an asymmetric outflow from the binary with velocities up to ±1200 km s~^. 
The origin of the outflow would be the emission source on the leading side of the 
accretion disc. 

Key words: Physical data and processes: accretion, accretion discs - stars: novae, 
cataclysmic variables - individual: RR Pic. 



1 INTRODUCTION 

Cataclysmic Variables (CVs) are close, interacting binary 
systems, comprising a white dwarf receiving mass from a 
Roche-lobe-filling late-type star. In absence of strong mag- 
netic fields, the mass transfer takes place via an accretion 
disc. RR Pic belongs to the subclass of CVs known as clas- 
sical novae, which have displayed a thermon uclear runaway 
outburst. The system has been discovered bv lSoencer Jonea 
il93ll) at maximum light in 1925 and, although it was of slow 
type, it is supposed to be in its quie scence stat e by now. 
The orbital period of 0.145025 days (IVoldll975^ places it 
just above the period gap of CVs. Vogt found the lightcurve 
dominated by a very broad hump, often interrupted by 
superimposed minima. He explained this behaviour by an 
exte nded hot spot region w ith an inhomogeneous struc- 
ture. iHaefner fc Met j il982l) however, explained their own 
time-resolved photometric and polarimetric observations to- 
gether with r adial velocity variations of the He II (4686 A) 
emission line JWyckoff fc Wehingeil[l977ll by suggesting the 
presence of an additional source of ra diation opposit e the 
hot spot. From high speed photometry. IWarnerl lll986ll con- 
cluded that during the 1970s (about 50 yr after outburst) 
structural changes have taken place in the system, result- 
ing in a more isotropic distribution of the emitted radia- 
tion. In addition, he has found evidence for a shallow, ir- 
regu lar eclipse, showing RR Pi c to be a high inclination sys- 
tem. iHaefner fc BetzenbichleJ 11199111 confirmed the general 
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change in the lightcurve of RR Pic as well as the presence of 
the eclipse, which they found to be very stable over several 
orbits. However, they point out that the additional emission 
source on the trailing side is still necessary to understand 
the phase relation between the radial velocity curve and the 
photometric lightcurve. 

Apart from this, the question arose whether RRPic 
is an in t ermed iate polar. Additional to the orbital period, 
iKubiakI lll984f) found a periodic modulation in the optical 
with a 15min period. He interpreted this as the rotation of 
the white d warf and concluded that RR Pic is an interme- 
diate polar. IHaefner fc Schoemb j il985ft however, repeated 
the high time-resolved photometry on a longer time-scale 
and could not find any sign of this short period. Since no 
15min period variation is found in X-ray measurements 
jBecker fc Marshal IT98lh either, they concluded that Ku- 
biak's variation was more likely a transient event in the disc 
rather than a rotating white dwarf. Also Warner's high- 
speed photometry does not reveal any period other than 
the orbital one. Hence no evidence remains that RRPic is 
an intermediate polar. 

With a quiescence magnitude mv = 12 mag, RRPic is 
a comparatively bright CV. Although it has been well ob- 
served photometrically, spectroscopic studies of the binary 
are sparse. In particular, no systematic investigation of the 
line profile has been undertaken so far, which could provide 
information on the accretion process in this system. In this 
paper we present new time-resolved spectroscopic data on 
RRPic and compare it with photometric data from the lit- 
erature. We establish a Doppler map of the Ha emission, 
study the isolated emission sources in the disc, and discuss 
the models of RR Pic's accretion disc. Finally, we analyse 
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Figure 1. The 19 individual spectra have been added up and di- 
vided by the mean continuum to produce this average, normalised 
spectrum. Only the part around Ho (6563 A) and He I (6678 A) 
is plotted. 



the high velocity wings of the Ha line and discuss possible 
interpretations. 



2 DATA 

The time-resolved optical spectroscopy was obtained on 
2002 February 23 with the Boiler & Chivens at the 1.52m 
ESO telescope on La Silla. Grating ^^20 together with a 1.5 
arcsec slit has been used to cover the wavelength range from 
5800 to 7500 A with a spectral resolution of 2.45 A FWHM. 
19 spectra were obtained, each with an exposure time of 
600 s. In total, we cover 3.8 h, which is slightly more than a 
complete orbit. IRAF has been used for the basic data re- 
duction including BIAS subtraction, flat-fielding, and wave- 
length calibration. No flux calibration has been performed. 
Two emission lines (Ha and He I) are present in the observed 
part of the spectrum. Both lines show a double peak profile 
in the average spectrum of all exposures (Fig. and strong 
variability in time (Fig. |5J . Hq has been used for further 
analysis in order to map the emission distribution in the 
system, thus gaining insight in the accretion structure and 
additional components. For this we used the tomography 
code by[Spruit (1998), wi th a MIDAS interfa ce replacing 
the original IDL routines (" Tappert et al.ll2002^ . 

Note, that this technique makes use of the variability 
of the line due to the binary rotation. We also expect a 
contribution to the line of the extended remnant shell of 
RR Pic whic h expands with a velocity around 400kms~l 
llRosino et alJrL982l) . However, since the shell is a circumbi- 
nary object, it is not effected by the binary rotation and 
any contribution to the line profile is a stable one. It does 
therefore not effect the here presented study. 



3 DISCUSSION 

3.1 Orbital period and zero phase 



The radial velocities have been determin ed by applyin g 
double Gaussian method as described in IShaften lll983 
evaluate the orbital period, two broad Gaussi ans have been 
fitted to each Ha line. The Scargle algorithm JScarglelll982l) 
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Figure 2. The time-resolved sequence of the two emission lines 
Ha (6563 A) and He I (6678 A) is plotted. Note that the scaling 
of the two plots is different to ensure maximal visibility of both 
lines. The little feature on the red side of Ha at epoch .714568 is 
an artefact due to a not properly subtracted cosmic ray. 



as implemented in MIDAS has been used to find the pe- 
riod in the radial velocities. As expected for data that cov- 
ers one orbit only, the periodogram shows a broad maxi- 
mum at P = 0.151 d. This maximum also i ncludes the value 
of P = 0.1450255(2) d, which IVogtl ^97^ has found from 
time-resolved photometry. We hence use this much more 
precise value for the further analysis. 

To compare the observed spectral features with the pho- 
tometric lightcurves in the literature, it is essential to syn- 
chronise the respective phases. Unfortunately, Vogt's mea- 
surements are too far away in time from our observations 
for the accuracy of the p e riod t o yield an unambiguous zero 
phase. However, iKubiakI (^8^ later showed the constancy 
of the period and ev en im proved the value slightly, in ad - 
dition IWarneJ fl986') and iHaefner fc Betzenbichieil lll99ll) 
agree on the stability of the eclipse. Bringing these infor- 
mation together, we finally made use of Wa rner's lightcurve 
S3464 from December 1984 J Warned 1198^ . determined its 
HJD of the eclipse as zero point, and added Kubiak's period 
to create the ephemeris 



eclipse = HJD2450064.441(2) + 0.14502545(7) • E, 



(1) 



with E being the number of cycles. These values are of suf- 
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Figure 3. Equivalent widtli of Ha plotted against photometric 
eclipse phase. Due to its breadth, the minimum at 4>e = 1 is 
unlikely to represent an eclipse. 



ficiently high accuracy to be extrapolated and compared to 
our observations. Defining the zero phase i^e = as the 
time of the eclipse we derive the orbital phase i^e for each 
observational point as 



HJDobs - 2450064.441 
0.14502545 



+ E 



(2) 



Our observation starts E = 15612 cycles after Warner's 
S3464 lightcurve. For the zero phase we hence derive an 
uncertainty of ao = 0.003 d or 0.02 phases. For comparison 
with Fig. ill the first observation (HJD = 2452328.578335) 
was obtained at phase (j>E = 7 • 10""", the typical time be- 
tween two exposures is equivalent to 0.053 phase units. 



3.2 Equivalent width 

In Fig. 121 the equivalent width of the Ha-line has been plot- 
ted against the orbital phase as derived from Equation [51 
A clear sinusoidal variation with a broad minimum around 
1.0 and maximum around 0.5 orbits is detected. This be- 
ha viour agrees with measurements of the equivalent width 
bv lHaefner ~B ctzcnbichlcr (199j). 

A compar ison w ith publis hed lig htcurves 
llHaefner fc Metd lll982ll. maxne^ l)l98(i) . and 
iHaefner fc Betzenbichlen (Il99lll l shows that the broad 
maximum in the lightcurves, which is partly eclipsed, does 
not coincide with maximal Ha emission, but rather with the 
broad minimum of the equivalent width. This indicates that 
it is probably caused by a region with enhanced continuum 
radiation, which is hence more optically thick, and therefore 
suppresses the Hq emission from this part of the disc. The 
maximal Ha flux instead comes from a region situated on 
the opposite side of the disc and coincides with the start of 
the minimum at cf>E = 0.7 (named No 1 by Warner) of the 
lightcurve. 

Important for the later computation of the Doppler 
maps is the fact that, within the errors of the equivalent 
width, no eclipse feature is found in the line. Hence, all data 
can be used for the Doppler tomography. Still, the varia- 
tion of equivalent width requires a careful handling, which 
is discussed in more detail in the corresponding section. 
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Figure 4. Diagnostic diagram for RRPic showing the param- 
eters of the radial velocity fit as function of the Gaussian sepa- 
ration. The phases of the radial velocities were calculated with 
respect to the eclipse ephemeris 4>e- The <po plot therefore gives 
the difference between <f>E and the zero point of the radial veloc- 
ity fit. The dotted line indicates the adopted correction resulting 
from the diagnostic diagram. 



3.3 Diagnostic Diagram 

While broad Gaussians have been sufficient to measure the 
orbital period, greater care must be taken to determine those 
radial velocities, which shall be used to study the line pro- 
file in more detail and to derive orbital parameters. As the 
motion of the white dwarf should be best represented by the 
extreme high-velocity line wings, which originate from the 
innermost part of the disc, one aims for a maximum sepa- 
ration of the two Gaussians without being contaminated by 
the continuu m noise. This is usually achieved by a diagnos- 
tic diagram (iShaf teJll98,'^ ) as shown in Fig. where the 
parameters of the radial velocity fit. 



Vr{<j)) ~ y ^ Ki sin 2-n(f), 



(3) 



are plotted as a function of the separation of the two Gaus- 
sians. The FWHM of each Gaussian was chosen to be 2A. 
The amplitude K\ shows an extremely steep slope to higher 
velocities with separation d and reaches a maximum value 
of about 350 kms~^. Although high values of Ki might be 
expected from a high inclination system such as RRPic, the 
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Figure 5. The radial velocities determined with the separation 
d = 34 of the double Gaussian are plotted against the orbital 
phase. The line indicates the best sinusoidal fit to the data. The 
semi— amplitude for this separation is Ki = 169.7it7kms~^. The 
offset in radial velocities can be explained by the system velocity 
7, the offset between zero phase of radial velocity and photometric 
eclipse is discussed in the text. 



value derived in this case is very improbable if not impos- 
sible. With Ki = 350 kms"^ and P = 0.145 d the distance 
for the white dwarf to the centre of mass is 0.9 solar radii. 
The maximum possible distance between white dwarf and 
secondary for a CV with this period is roughly 1.6 solar 
radii. Hence, the centre of mass would be situated nearer 
to the secondary than to the primary, which would imply 
a mass ratio M2/M1 > 1, and does not allow for a stable 
mass transfer. Additionally, we point out that the maximum 
values of Ki are reached at very high separations, which are 
already far away from the line centre in a region of very 
low S/N. We conclude that even at high separation, the line 
is still affected by additional emission sources, which then 
must be large and reach either deep into the central part of 
the disc or do not belong to the disc at all, as discussed in 
section 13.51 

As a consequence of the steep increase of Ki, the or- 
dinarily plotted ratio (j{Ki)/K\ stays nearly constant and 
can not be used to indicate maximal possible separation. We 
instead plot the error of the sinusoidal radial velocity fit a to 
give an idea of the uncertainties. The strong variation that 
starts around d ~ 34 might indicate the point of optimal 
separation and favours a semi-amplitude Ki ~ 170kms~^. 
In Fig. |S] the radial velocities for d = 34 are plotted against 
the orbital phase. 

After t he initial disto rtion due to the double peaked 
line profile (iTappertjllOOgl) the zero phase settles at a value 
around 0.05 between separations d 10 and 20. A sim- 
ilar value has also been found by iHaeftier &: Met j lll982ll 
who co mbined their lightcurves wit h radial velocities mea- 
sured bv lWvckoff fc Wehingeil (^T^) . This value for the zero 
phase is, however, very likely to be infiuenced by additional 
emission sources, which give a strong bias at low separa- 
tions. Wyckhoff & Wehinger state that they derived the ra- 
dial velocities with "classical methods", and therefore very 
probably used the total line for their measurements. Hence, 
their zero phase can be assumed to be biased in the same 
way as our values at small d. 



Emission 
from 

leading side 





Emission 
from 

trailing side 



Secondary 




Eclipse at 
4> <0 



4> = 



Figure 6. The sketch visualises the geometry of an eclipse (left 
side) before the superior conjunction of the white dwarf {4'r = 0, 
right side). The cross marks the centre of rotation. The eclipsed 
part has to be an emission source on the leading side of the disc. 



For higher separations (d ~ 30 and above), which are 
usually believed to more likely refiect the motion of the white 
dwarf, the zero phase reaches 0o = 0.16 and stays almost 
constant afterwards. It therefore very likely corresponds to 
the superior conjunction of the white dwarf and gives an 
unambiguous input value for the subsequent Doppler map- 
ping. Note that, since the orbital phase was selected in such 
a way that = for the eclipse (see Equation the zero 
phase 00 indicates the difference between the phase of the 
eclipse and the zero phase of radial velocities; i.e. the eclipse 
appears 0.16 orbits before the superior conjunction of the 
white dwarf (see also Fig. (SJ. This can only be explained 
with the assumption that not the disc as a whole is eclipsed, 
but rather an emission region on the leading side of the disc 
(Fig.li. 



3.4 Doppler map 

For the interpretation of the phase dependent line pro- 
files, the Doppler tomography method as developed by 
iMarsh jTli ornc ( 1988) has been applied. A Doppler map 



displays the fiux emitted by gas moving with the 



velocity {vx,Vy). This deprojection is achieved due to the ro- 
tation of the binary system, as the projection of the velocity 



follows a sinusoidal radial velocity curve. Hence, the 



line profiles can be transformed into a map I{vx,Vy) without 
any assumptions about the actual velocity field. 

To correct for the variation of the equivalent width, the 
input spectra have been normalised by the Ha emission line 
fiux. The intensity values in the Doppler map are hence to 
be interpreted as relative fiux values only. Furthermore, we 
corrected the orbital phase of the input spectra by the zero 
value 00 = 0.16. 

The final Doppler map is given in Fig.|7|and shows the 
distribution of Ha emission sources in velocity coordinates. 
Apart from the disc itself, which is clearly visible, two signif- 
icant additional emission sources can be seen. The orienta- 
tion of the map is chosen in such a way that the phase angle 
(l)r with respect to the zero point of radial velocities, is zero 
towards the top and increases clockwise. The smaller of the 
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Figure 7. Doppler map of RRPic showing the distribution of 
HcK-sources in velocity coordinates. The cross marks the centre of 
rotation, which corresponds to the centre of gravity of the system. 
The map is oriented in such a way that the phase angle referring 
to the zero point of radial velocities is zero towards the top and 
increases clockwise. The contour lines refer to values 75, 85, and 
95 and indicate the relative brightness of the emission sources. 
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Figure 8. V/R and its residuals is plotted against the orbital 
phase referring to radial velocity zero point. Error bars have been 
determined by a Monte Carlo simulation. Two orbits are plotted, 
the original data (o) with error bars between phases 0.0 and 1.0, 
and the reconstructed (o) one with error bars between phases 1.0 
and 2.0. The latter points have additionally been connected via 
dashed lines. 



two emission sources is situated at a phase angle 4>r ~ 0.8, 
a place that is usually identified with emission from the hot 
spot. The stronger emission, however, arises on the opposite 
side, 180° away, at (j>r ~ 0.3. 

The quality of the Doppler map can be judged in sev- 
eral ways. In this paper a quantitative approach is made to 
compare the shape of the emission line in the original and re- 
constructed spectrum. The asymmetry of the line profile can 
be measured by a V/R plot. Here, the line profile is divided 
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Figure 9. The observed Ha emission line profile of RRPic has 
been plotted against the phase. The spectra have been flux nor- 
malised as described in the text and the wavelength transformed 
into radial velocities with respect to Aq = 6562.8 A. Intensities 
are logarithmically colour coded to favour the display of the faint 
outer wings. The little enhancement at phase 0.83 and 900 km s~^ 
is due to a not properly subtracted cosmic ray. 

in two halves, and V/R defined as logarithm of the ratio of 
the fl uxes in the blue (violet) and red half jTappert et alJ 
l2002h . In Fig.|H| V/R of the original and reconstructed line 
profile is plotted against the orbital phase. The residuals 
given below show no sign of periodicity, hence indicating 
the good quality of the fit. 

3.5 High velocity wings 

In Fig. 1^ the changing line profiles of the flux normalised 
Hq emission line are displayed as an image with emphasis 
on the faint outer parts of the line. One can clearly see that 
the line wings reach up to velocities of 1200 km s~^ and vary 
with the orbital phase. 

The classical explanation for this feature is that it is 
caused by gas circling the white dwarf at velocities of about 
1200 km s~^. However, there are some problems with this 
picture. According to Fig. I^l the gas not only has these high 
velocities, but also changes its velocity by the same mag- 
nitude with phase; it is completely on the side of negative 
velocities at phase (pr ~ 0.25 and on the side of positive ve- 
locities at <l)r — 0.75. As the high velocity gas is supposed 
to be in the inner part of the disc, this indicates a semi- 
amplitude Ki « 1000 kms~^, which is a very unlikely value, 
especially as the distance of the white dwarf from the centre 
of mass in the Doppler map indicates a much smaller Ki . 

Similar high velocity wings have been found in the 
Balmer lines of t he SW Sex-type systems, e.g. LSPeg 
jTavlor et al.lll999ri . They explain the wings and their vari- 
ation with a mo dified disc-overflow model similar to the 
one developed bv lHellier fc RobinsonI il994) . There are cer- 



6 L. Schmidtobreick, C. Tappert, I. Saviane 



tainly similarities between RRPic and SW Sex-type stars, 
as the orbital period is between 31i and 4h and the radial 
velocity phase lags behind the eclipse phase (which we have 
explained by an eclipse of the emission source on the lead- 
ing side). However, the dominant absorption feature around 
phase 0.5, which is typical for SW Sex-type stars and prob- 
ably related to the disc-overflow is not present in RR Pic. 
Also, if all accreted material were deflected over the disc, 
no hot spot could be present. However, some SW Sex-type 
stars are known to show emission sources on the trailing side 
of their discs. Hence the possibility of part of the accretion 
stream flowing over the disc can not be ruled out for RRPic. 

A different explanation is that the high velocity gas does 
not stay inside the system but is actually ejected in the form 
of an asymmetric outflow, with maximum projected velocity 
of 1200kms~'^. Taking into account the proposed inclination 
i ~ 65° (Jia efne r fc Mctz 1982) the real outflow velocity has 
to be around 1350 km s^^. This would explain the varia- 
tion of the high velocity component as a lighthouse effect: 
the outflow emerges from one part of the disc only and is 
swept around with the rotation of the system. In this pic- 
ture, both the high velocity itself and its varying amplitude 
are explained in a consistent way. Including the information 
from phase (j)r = 0.25 of maximal negative velocity (outflow 
towards us), the origin of the outflow coincides with the po- 
sition of the isolated emission source on the leading side of 
the accretion disc. See Fig. llUl for an explanatory sketch. 

Additional evidence for the presence of an outflow in 
RRPic comes from the mapping of its remnant shell, which 
shows a bipolar structure in the form of an equatorial ring 
with some tails orthogonal to it jGifl fc 0'Brienlll998l) . Al- 
though this map shows the shell far outside the binary, the 
ring might well be fed by the asymmetric outflow from the 
binary. The fact that no P Cygni profiles hav e been found 
in lUE spectra of RRPic jRosino et al.lfl9 82|) has been in- 
terpreted as RRPic having no outflow and that the shell is 
already decoupled from the binary. The asymmetric outflow 
can be considered as optically thin as it is visible in emis- 
sion. Any P Cygni absorption has hence to be weak, and 
can furthermore only be produced during the phase when 
the outflow is directed towards the observer. It would hence 
smear out over normal exposure times, and is therefore not 
in contradiction with the observations. 



4 CONCLUSION 

We have shown via a Doppler map that two isolated emission 
sources are present in the disc of RRPic. Combining the 
radial velocities with published lightcurves we come to a 
rather uncommon interpretation of the photometric eclipse 
as an eclipse of the leading side emission source (see Fig. (SJl 
rat her than an ec lipse of the hot spot as suggested previously 
by|Kubia3(Il983). 

Putting all th e information together, we conflrm the 
working model of iHaefner fc Met j Jl982h for RRPic in 
preference to the one of IVoetl Jl975l) . Although, as dis- 
cussed in section 2] the superior conjunction of the white 
dwarf takes place 0.1 earlier in phase, this does not change 
their basic interpretation. Warner's (1986) doubts of the 
model were based on the fact that no isolated emission 
from the leading side had been found in other CVs. Mean- 
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Figure 10. The sketch visualises our model of the emission 
sources in RRPic and the outflow coming from the emission re- 
gion on the leading side of the disc. During a complete orbit of 
the binary rotation, the projected velocity of the outflow changes 
between negative (flow towards the observer) and positive (flow 
away from observer) velocities. 

whi le many systems with this structure have been found 
fsee lTappert fc Hanuschi^ ll200 for an overview), and this 
argument is no longer valid. 

Although common, no convincing explanation for iso- 
lated emission sources on the leading side of the disc has so 
far been presented. Among other possibilities, the presence 
of spiral shock waves has been discussed as possibl e expla- 
nation for similar features in other CVs (see e.g. ISteeghs) 
i200ll) l. Due to their still enhanced mass transfer, long- 
period old novae should possess extended hot discs, and 
are therefore indeed candidates for possessing spiral shocks 
jBoflinll200l . However, although the emission sources in 
the disc of RR Pic have a clearly elongated structure, nei- 
ther time nor spectral resolution of our data is high enough 
to confirm this picture for RRPic. 

The analysis of the high velocity wings of the Ha line 
yields two possible explanations for these wings and their 
variability. They may originate from a partial disc-overflow, 
or from an asymmetric outflow from the accretion disc, i.e. 
from the emission source on the leading side of the disc. 
In this connection between the outflow of the binary 

and the ring-like structure of RR Pic's shell seems plausible. 
Deep, high resolution Ha images as well as high resolution 
spectroscopy of some metal lines might help to clarify this 
picture. 
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